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_ -, ' The deformation of a nucleon embedded in various finite nuclei is considered 

■ by taking into account the distortion of the chiral profile functions under 

the action of an external field representing the nuclear density. The baryon 
. charge distribution of the nucleon inside light, medium-heavy and heavy 

nuclei is discussed. The mass of the nucleon decreases as it is placed deeper 
| inside the nucleus and reaches its minimum at the center of the nucleus. We 

' discuss the quantization of non-spherical solitons and its consequences for 

^ ■ the mass splitting of the A states. We show that bound nucleons acquire an 

i . intrinsic quadrupole moment due to the deformation effects. These effects 

are maximal for densities of nuclei about p(R) ~ 0.3 4- 0.35p(0). We also 
point out that scale changes of the electromagnetic radii can not simply be 
described by an overall swelling factor. 
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I. INTRODUCTION 



The possible modification of nucleoli properties in the nuclear medium is currently a 
much discussed topic in low energy hadron physics. One way to consider such problems is 
to describe the nucleon as a topological soliton and then study the influence of baryonic 
matter on the properties of such solitons, making use of the independent particle picture 
which has been so successful in describing many properties of nuclei. There exist already 
some works where nucleons described as Skyrme-type solitons embedded in infinite nu- 
clear matter have been considered 0j2[].f3 The results of these studies were in qualitative 
agreement with experimental indications and with results of other authors using different 
approaches. On the quantitative level, however, there is a too large renormalization of the 
nucleons' effective mass in nuclear matter. The quantitative value of this renormalization 
is about ~ 40% for the normal nuclear matter density, po = 0.5m,,: 3 [3J. It is therefore dif- 
ficult to relate this modification of the nucleon self-energy in the medium to the nucleon 
mass in free space because such calculations consider only the special case when a nucleon 
is placed in the center of the heavy nucleus. The infinite nuclear matter approach allows 
one only to consider properties of nucleons placed near the center of heavy nuclei where 
the density is constant. One can expect that taking into account non-spherical effects, i.e. 
deformation of the skyrmion in the finite nucleus, would improve the results also on the 
quantitative level. Density changes play an important role when the nucleon is placed at 
sufficiently large distances from the center of nucleus. 

The role of deformation effects (for an early work see ||) has already been investi- 
gated in connection with the nucleon-nucleon interaction in ref. @. The results presented 
in this work demonstrated the necessity of non-spherical deformed solutions of the field 
equations to obtain a reliable description of the NN potential. A similar approach to the 
deformation for two interacting skyrmions has been used in the work 0, where changes 
in the nucleon shape are investigated by letting the nucleon deform under the strong 
interactions with another nucleon. On the other hand, rotational and vibrational excita- 
tions of deformed skyrmions have also been considered for baryon number equal to one 
(B = 1) system ||. In a series of papers [pHTTH axially symmetric few-skyrmion systems 
with baryon number B larger than one were investigated, because spherically symmetric 
configurations in the skyrme model do not give bound skyrmionic systems Jl2| even for 
the solutions of the model within the finite nuclei as considered in ref. jTB| (for similar or 



related work, see also ref. |l4p. In the works the deformation is introduced by 



a modification of the unit spherical vector and additionally through the modification of 
the profile function || in the expression of the ansatz for the chiral field (for more precise 
definitions, see below). The results of all investigations for B > 1 systems lead to the 
conclusion that the system has to be deformed to obtain the ground state. Some other 
papers that deal with deformation effects and non-spherical configurations are collected 



in refs. [TT5 



*For an earlier and somewhat different approach to density effects in Skyrme-type models, 
which leads to similar results, see ref. ||. Some other works concerned with nuclear matter 
aspects of the Skyrme model are collected in ref. 
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In the present work we consider properties of the deformed nucleon embedded in light, 
medium-heavy and heavy nuclei. Deformation effects are introduced by the distortion 
of the profile function of chiral field under the action of the external field (which pa- 
rameterizes the baryonic density within a given nucleus). In contrast to most previous 
investigations, we do not consider a constant nuclear density as appropriate for nuclear 
matter but rather the distance-dependent density as given for finite nuclei. We calcu- 
late the modifications of the nucleon properties in finite nuclei, and their dependence on 
the distance between the topological center of the skyrmion inside the nucleus and the 
geometrical center of the nucleus under consideration will be considered. 

The paper is organized in the following way. After the formulation of the problem 
in section |J, we discuss the classical soliton mass and distribution of baryon charge for 
a skyrmion inside finite nuclei in section JTJ. Based upon the results for the mass and 
baryon charge distribution of the deformed classical soliton, we present in section |V| 
the quantization of axially symmetric systems with baryon number B = 1. Section [V] 
is devoted to the quantum properties of the solitons, in particular the effective mass of 
nucleons in finite nuclei and the consequences for A mass splittings. We also consider 
the density dependence of the nucleons electromagnetic radii and calculate the intrinsic 



quadrupole moment of the deformed nucleon. Finally, in section [VI] the conclusions and 
an outlook for further investigations are given. 



II. FORMULATION OF THE PROBLEM 
A. Lagrangian and ansatz 



Our starting point is the medium-modified Skyrme lagrangian 

-l-(f)(f)-i^(w)(w t) 



1 r Tr[L M , Uf + S^La s (f)Tr[U + - 2] 



(1) 



32e 2 16 



where L M = U + d fl U, F n = 186 MeV is the pion decay constant, e is a dimensionless 
parameter and m n = 138 MeV is the pion mass. The medium functionals a s (r) and a p (r) 
are expressed via the functionals Xs(r) and Xpif) °f s-wave and p-wave pion-nucleon 
scattering lengths and the nuclear density p(f): 

a p (r) = 1 - X p(r) , a s (f) = 1 + ^ . (2) 

Here, Xp{^) nas t ne meaning of the pion dipole susceptibility of the medium. The influence 
of the medium enters through these functionals (for explicit expressions, see section [HC] ) -H 



^Such an approach has also been considered in next-to-leading order chiral perturbation theory 
in which the pions are chirally coupled to matter fields, like e.g. nucleons, see Ref. fllqj. 
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Note that we employ here the same approximations as in |L]. Since the fourth order term 
can be obtained from the infinite mass limit of p meson exchange, a possible density 
dependence of that term could enter via in-medium modified p meson properties as it 
was done e.g. in |||. Similarly, one could add a sixth order term ~ B^B^ arising from 
uj meson exchange or an symmetric fourth order term due to scalar mesons. We do not 
consider such terms here. Also, since the empirical information on possible changes of 
vector meson properties in the medium has not yet lead to an established picture, we 
treat the coupling e as a density independent number. This could be released in more 
sophisticated model but goes beyond the scope of the present investigation. For a more 
detailed discussion about how the specific form of the Lagrangian in eq. (0) arises, we 
refer to ref. Q. 

The ansatz for the chiral field U has the hedgehog form 

(3) 

for the skyrmion placed at the center of the nucleus. 

When the skyrmion is placed at a distance R from the center of the spherical nucleus 
(see Fig. [I]), the chiral field U can be written in the form 

U(f— R) = exp 

where N is the unit isotopic vector correlated with the unit space vector (r — R)/\r — R\. 
In the more general case, N is a functional of some angular functions 0(9), tp) PH|JTT[] . 
The dependence N = N[Q(9), <&(9, p)} and the dependence of the profile function F on the 
direction of the radial co-ordinate space vector r — R characterizes the deformation of the 
skyrmion. As it was mentioned above, in the works fl-pl such deformation was taken into 
account when considering properties of multi-skyrmionic systems. These modifications are 
dictated by the topological properties of the model, in particular the $(6*, p) dependence 
is dictated by the constraint that the baryon number is larger than one, B > 1. Here we 
consider the B = 1 case and simply set $(9, p>) = p. Consequently, we can represent the 
vector N in the following form 

N = {sin6 (9) cos p, sin9 (9) sin p, cos @(9)} . (5) 

However, this is not the only deformation effect. We also modify the profile function F 
by allowing for a 9 dependence, 

F — F(\f— R\, 9) . (6) 

Note that any azimuthal <p dependence of F is removed due to the axial symmetry of 
the system (see the next subsection). Obviously these deformation effects also exist in 
the case when a skyrmion is placed at the center of nucleus, R = 0, but in this case the 
deformation has only a radial dependence. In the previous work |]J this dependence was 
determined by solving the equations for the profile function of the ansatz (|3]) in nuclear 



U{r) 



exp 



it ■ r 



F(\r\) 



if-N\ I F{r-R) 



(4) 
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matter, F(r) — > F*(r). It is clear that in the present case the deformation has in general 
both a radial and an angular dependence. The angular dependence will come into play 
when the skyrmion is located at some distance from the center of the nucleus. The radial 
part of the deformation represents the breathing mode of the skyrmion and the angular 
part represents its deviations from the spherical form. Finally, taking into account our 
choice of N and F (eqs. @ and @), the following ansatz for the chiral field U arises 

U(r-R)=exp[iT-N(&(6),<p)F(\r-R\,6)] . (7) 

This form will be used in the present work. As a result some modifications beyond the one 
induced by the breathing mode will appear. For example, one expects that an intrinsic 
quadrupole moment of the skyrmion emerges. This will be discussed in more detail below. 

Before discussing the in-medium skyrmion properties, we point out that throughout 
this paper we neglect the effects due to Fermi motion of the solitons (nucleons). In a 
strict large-iV c counting, such effects are suppressed. However, since we do not expect 
deformation effects to be large, one can doubt the validity of this formal suppression. At 
present, it is simply not known how to include such effects systematically because one 
not only would we have to boost the solitons (where the boost velocity could be obtained 
in the local density approximation) but also the coupling to the rotational modes would 
considerably complicate the quantization procedure to obtain physical states. While the 
neglect of Fermi motion might to some extent alter the quantitative results of our study, 
we believe, however, that the qualitative ones would remain unaltered. 



B. Mass functional in the static case 



The mass functional in the static case can be obtained from the lagrangian ([!]) by 
using the ansatz (|7|) 

M(R) = -J £ static [C/(f- R),f]d 3 f. (8) 

To simplify the calculations one introduces the variable r' — r — R putting the origin of 
the coordinate system O at the topological center of the skyrmion O' (see Fig. |IJ). Then 
the expression for the mass functional can be written as 

M{R) = - J C static [U(f'),f',R]d 3 f'. (9) 

We note that the distance \R\ between the skyrmion and the nucleus is a matter of 
convention. It depends on how one determines the center of the soliton. If the skyrmion 
is deformed its mass and topological centers may not coincide. In this paper we determine 
\R\ as the distance between the center or the nucleus under consideration and topological 
center of the skyrmion under consideration. 

From figure [l] it is easy to see that the system has an axial symmetry. One can choose 
the direction of the z' axis along the direction of the vector R so that the integrand 
in eq. (§) has no azimuthal ip' dependence. Consequently the profile function F only 
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depends on the variables r' and 9'. Therefore, the lagrangian (JJ) in the static case can be 
written in the form 

Static = —[{diNj ■ diNj) sin 2 F + diFd^a^r', 6\ R) 

o 



1 

1 

2^2 



[{diN ■ diN) 2 - (diN ■ d 3 N) 2 } sin 4 F 

[{diN ■ diN)djFdjF - (diN ■ djN^FdjF] sin 2 F 



m 2 F 2 

— ^(1 - cos F)a s (r',e',R) 



(10) 



which (after some algebraic manipulations) leads to following mass functional^ 

, wx,s t^F-e 7 J \ sin 2 6 sin 4 F 9 1 9 ^ / sin 2 6 ^ 9 \ _ „ ~, 

e j j sin v r A 4 \ sin v I 

oo^ x ' 



+ 



sin 2 sin 2 F a p (f, 9, R) 
sin 2 Of 2 4 



(ii) 



+ 



V sin' 



4 



F? \sm 6d9dr. 



dF 



OF 



Here F? = -r— -, = — — , = and we have introduced the dimensionless variables 
or oO 

f = eF n r\ R = eF^R. 

The corresponding baryon number of the skyrmion is given by 



OO 7T 

B= / / F f & 6 sin 2 F sin QdOdf . 

o o 



(12) 



The static soliton solution is found by minimizing the functional (11). For doing that, 
one has to solve a system of coupled differential equations 



f(Fff, F e e, F?, F e , Qg, F, 0) = 
g(Qee,Qe,F f ,Fe,e,F) = , 



(13) 



for the functions F(f,9) and 0(0) subject to the boundary conditions determined from 
the baryon number condition B = 1 (|i~2|). To simplify the numerical computations we 



* Hereafter, we drop the prime on the variable 6 for convenience. 
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perform various approximations. First, consider the profile function F(r) for an unde- 
formed (spherical) skyrmion. It is well known that F(r) can be parameterized to a very 
good approximation as 



2arctanl^-l (14) 



for the free skyrmions with baryon number B — 1. The value of the parameter r$ is 
obtained from the constraint of having the minimum mass functional for the free skyrmion. 
This parameterization in principle can be used for a bound skyrmion and we apply this 
ansatz in constructing our profile function (despite some short-coming, see section |III| ). 
Also, from previous calculations we know that the dependence of the profile function on 
the density is weak, see e.g. [Q. Renormalization of the mass functional (|ll|) mainly occurs 
via the medium functionals a v and a s . Furthermore, from Fig. ([I]) we see that changes of 
9' mainly affect the value of medium functional^] and following the results from ref. 0, 
we assume changes of profile function to be small. Using these facts we can represent our 
profile function as 



r 



F = 2 arctan j ^ [1 + ^ cos 9 + 72 cos 2 9 + 73 cos d 9 + . . .] j , (15) 

where the cos functions are chosen to maintain periodicity in 9. For the function we 
choose the parameterization 

= 9 + 6 1 sin 29 + S 2 sin 49 + S 3 sin 69 + . . . (16) 

in order to avoid singularities from factors like sin / sin 9 in the expression of the mass 
functional ([□]) ||T0| , pT| . In eqs. (|15|) and (|16"1) , rs, 71, 72, 73, ••• and Si, 62, S3, ... are 



variational (or in other words deformation) parameters. 

Our aim is the minimization of mass the functional ([11]) using the expressions (|T5|) 
and fllBT) for the profile function. Before doing that, we have to specify the medium 
dependence that enters through the s- and p-wave pion-nucleon functionals Xs and x v , 
respectively. 

C. Medium functionals and parameterization of the density 

The medium functionals \s and Xp introduced in Eqs. ([[]), (fj) have the form 
Xs(r,e,R) = -A7rr]b oP (r,9,R) , 

»' *> = T^TUT ' < f ' % = Am ' AfAk) , (17) 

l + g' K(r,9,R) V 



* We remind that a p and a s are functionals of the medium density p{r). 

**We consider here symmetric nuclei with p n = p p . For asymmetric nuclei, additional isovector 
effects will be present but are not considered here. 
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where rj = 1 + m^jM^ ~ 1.14 is a kinematical factor, = 938 Me V the mass of 
the nucleon and g' = 1/3 is the Migdal parameter which takes into account the short- 
range pair correlations of the dipole centers. The empirical parameters &o — — 0.024m" 1 , 
Cq = 0.21m~ 3 [17 can be taken from analyses of low energy pion-nucleus scattering data. 



We choose the following parameterization for the density of the finite spherical nuclei 
considered here 



p(r) 
p(r) 



A - I 



'A- I 



1 + 



Z - 2 / r 



3 Vr 




exp <^ — ? r > for 4 < A < 20 



for A > 20 . 



(18) 



In the actual calculations (minimization procedure) we have considered the nuclei 12 C, 
16 0, 40 Ca, 56 Fe, 198 Au and 208 Pb. The first two represent light, the next two medium- 
heavy and the last two heavy nuclei. The parameter ro depends on the type of nucleus 
considered. Its value is 1.31 fm for 12 C and 1.76 fm for 16 0. Z is the charge of the nucleus 
and po = 0.5m~ 3 is the normal nuclear matter density. Furthermore, R! = 1.2v4 1//3 fm 
and a = 0.58 fm. A is the baryon number of the nucleus and the factor (A — 1) j A allows 
one to consider properties of a given nucleon which is one of the nucleons of the nucleus. 
For the accuracy of our calculations, such a simple rescaling of the nuclear density should 
be sufficient. The density is pameterized as 



p(f,6,R) = p(^f 2 + 2rRcos6 + R 2 ) 
as dictated by the axial symmetry of the system. 



(19) 



III. MASS OF THE CLASSICAL SKYRMION AND THE BARYON CHARGE 

DISTRIBUTION 

For the input parameters of the Skyrme model we use F n = 108 MeV and e = 5.265 
in order to reproduce the free space masses of the nucleon and the A as in the seminal 
work ||19|| . One could also work with the empirical value of F n and adjust e to the iVA 
mass splitting. In that case, however, the nucleon mass would be too large, but relative 
changes would be affected considerably less 0. For definiteness, we prefer here to work 
with the physical nucleon mass. 

Our calculations have been performed using different sets of the variational parameters 
r s, 7i, 72 , 73, • • • , Si, 62, S3, ... . We represent the skyrmion size parameter r$ as 

rs = r f s - Sr s , (20) 

where is the value for the skyrmion in free space (r| = 0.6 fm). 

The minimization procedure of the mass functional (|11|) showed that considering the 
set of parameters 
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f 7i = 0, if i > 3 , 
\^ = 0, if j > 2 , 



(21) 



is sufficiently accurate. The strength of the parameter 7$ decreases with increasing index 
i. For example, 73 is one order of magnitude smaller than 7! at the distances where 
deformation effects are maximal. The parameter 5j decrease very fast with increasing 
index j. For example, 62 is one order of magnitude smaller than 5\ and £3 is almost zero. 
For controlling our results and keeping a high accuracy, we have performed calculations 
with the following set of parameters: 



For the reasons just mentioned, we only show results of 71, 72 and 8\. 

The dependence of the skyrmion mass from the distance between the centers of the 
nucleus considered and the soliton (which represents one of the nucleons within the nu- 
cleus) is shown in Figs. [2|, ^ |4], for the light, medium-heavy and heavy nuclei, in order. 
The bands labeled "1" in these figures represent the soliton mass normalized to its free 
space value and for convenience we display by the bands denoted "2" the nuclear density 
normalized to its value at the center of the nucleus. These bands are shown in the top 
panels of these figures. The behavior of the soliton mass is similar for all cases when the 
soliton is placed within various nuclei, i.e. the mass increases with falling medium density 
and approaches its free space value at the border of the nucleus. It has its minimum value 
at the center of the nuclei. This is even the case for the light nuclei, which do not have 
their maximum density exactly in their center. 

The deformation parameters which represent the deviations of the nucleon shape from 
the spherical one are presented in the bottom panels of Figs. 0, |3], f|. The bands denoted 
"3" represent the dependence of the relative value of the radial deformation parameter 
5rs/rs from soliton location inside nucleus and the bands "4" and "5" give the dependence 
for 71 and 72, respectively. The band "6" represents the parameter 8\ multiplied by factor 
10. The other deformation parameters are very small and we do not show them here. The 
leading parameter 71 has its maximum value at the distances where the density of the 



nuclei falls off by -20% for 12 C and 16 0, ~ 40% for 40 Ca, 56 Fe, 198 Au and 208 Pb, whereas 



the radial deformation parameter 5r$ has its maximal value at the center of the nuclei. 
The parameters 72 and S\ have maximum values where the density of the nucleus falls off 
~80%. In these regions, the density gradient is maximal leading to the largest angular 
deformation. From this one concludes that the skyrmions inside nuclei is in a deformed 



The topological baryon charge distribution for a skyrmion embedded in the nuclei 12 C 
and 56 Fe is shown in Figs. [5|, |6], respectively, where figures a) represent the baryon charge 
distribution in the y, z plane and figures b) give its projection onto this plane. Each area 
at the bottom figure represents the projection of surfaces of equal baryon charge density. 
The point y — z — indicates the topological center of the skyrmion and the center of 
nucleus is situated at the negative side of the z axis. One can see that opposite directions 
in the z direction are not equivalent in this plane. Stated differently, going for the same 
distance in opposite directions in z gives different values of the baryon charge density. 




(22) 



state. 
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We remark that the holes at the origin of the system (y = z = 0) are caused by the 
parameterization of the profile function (see eq. Q15D ). In this approximation, Ff — > 
as f — > 0. This derivative takes a finite value if one directly solves the equation coming 
from minimization of the mass functional (|TTD not using the ansatz eq. (|l4l), thus no such 
holes would appear. However, the effects discussed so far and in the following are genuine 
and independent of the (dis) appearance of these computational artifacts. The light area 
at positive values of the axis z is the area of maximal density. The mass distribution, 
i.e. integrand in the expression ( |TTD shows a behavior similar to the one of baryon charge 
distribution. 

An exotic behavior of the mass and baryon charge distributions is observed when the 
skyrmion is placed at sufficiently large distances from the center of the nucleus, where 
density changes are significant (see Figs. ^, One can see that the maximum of the 
baryon density is shifted away from the center of the nuclei. When the skyrmion is located 
at the center of nucleus, the baryon density distribution is symmetric and thus the shift 
becomes more effective with increasing of R. The shift of the density maximum from 
the topological center of the skyrmion in opposite direction to the direction of center of 
the nucleus can be explained in the following way. The mass of skyrmion decreases in 
the nuclear medium and if the density is decreasing, the skyrmion mass increases. While 
we are going from the center of the nucleus to its surface the density fall-off and mass 
concentration of the skyrmion occurs at the side of the smaller nuclear density as seen 
from its topological center. A similar behavior of the baryon charge density for other 
nuclei is observed, so we do not present it here. 

Finally, at the end of this section we remark that rotations of the figures b) around the 
z-axis give three dimensional surfaces with an equal baryon charge distribution. From 
this we conclude that the soliton has two equal components of the moment of inertia 
which are not equal to the third one, I x = I y ^ I z . We will use this fact in the next 
section. 



IV. QUANTIZATION PROCEDURE FOR NON-SPHERICAL SOLITONS 

As it is well known, baryons emerge in the Skyrme model as quantized solitons [TP 



This quantization is performed in terms of an adiabatic rotation of the soliton and gives 
rise to the appropriate quantum numbers for the emerging physical states. In such an 
approach, the grand spin K, i.e. the sum of angular momentum J and isospin T, is 
conserved. This leads to a tower of states with J = T. It is well known that in this case 
only the nucleon (J = T = 1/2) and the A (J = T = 3/2) have physical meaning. All 
other particles (higher rotational states) should be considered artifacts of the model and 
will be discarded. 

To quantize the solitons of our model we use the standard canonical quantization 
procedure performing time-dependent rotations in isospace and in co-ordinate space 

U(f") = a(f') + ir-$(f') — ► U(x,t) = a(x) + %T%j (t)<P(x) , 

(23) 
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where the matrix Ijj specifies the rotation in isospace and R^ the one in co-ordinate 
space ||. The parameters of the rotation matrices ly and Rjj serve as collective variables 
describing the rotational degrees of freedom with respect to a body-fixed system of refer- 
ence. The latter has its origin at the topological center O' of the skyrmion and points in 
the direction of the distance vector R from the center of the nucleus to the center of the 
skyrmion (see Fig. |l|). 

The quantization procedure used here is similar to the standard one Jl9| but now the 



axial symmetry of the system should be taken into account |q-|Tl|]. The angular velocities 
of isorotations u>i and space rotations fli can be determined from the respective equations 



(24) 



where dot means differentiation with respect to the time component, and is the totally 
antisymmetric tensor in three dimensions. 

Using now the ansatze (|7]), (^3j) in the lagrange density ([I]) and integrating over the 
whole space one gets following expression 

L(R) = J C(R)d 3 f' = -M(R) + H 2 + ^)i£?(fl) + (ilj + S%)lgZ\R) 

-2(0^ + uj 2 n 2 )I^\R) + (a* - n 3 ) 2 I^\R) , (25) 



7T 



2e 3 F. 



/ /sin 2 F (3 + cos20) ^ + ^(1 + 4F 2 ) 
o o L ^ ' 



+ sin 2 F 



'sin 2 © 
sin 2 9 



+ 6 2 cos 2 6 



sin OdOdr 



(26) 



7T 



2e 3 F w 



o o 



r 2 sin 2 cos 2 9 
T sin 2 

. sin 2 9 



9 „ (f 2 (3 + cos 29) sin 2 9 „ 
sin 2 F+ — + ^ 



4 



sin 2 # 



sin 2 F 



F 2 



+ 7T f2 + 4 — ^7T sin2 F sin(9 
9 1 sin^ ' 



j. 2 



,sin 2 6 



+j ( 2[2-cos2(6-0) + cos2^]^ I ^sin 2 F + f :2 sin 2 (e 



sin 2 9 



-9) F 



+i sin 2 F (r 2 + 2(3 + cos 20)^? sin 2 f\ Ol 



sin 2 9 



-— sin 2 F[-3 + cos 2(6 - 9) ]F?Q 2 g 



sin 9d9dr 



(27) 



OO 7T 





'F 2 f 2 , „ A cos 

f + ¥< 1+4 ^hw sm2e 



/ ^2 in 2 q _ r ^ _ 

\ 4 sm^ 9 



COS6 1 ^ 99 ^ ^ 

H -0^ sin F sin 20 

2 sin 



sin 9d9dr , 



(28) 
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' :; ; 1 ' " 1 - y y sin 2 F I T — + 2F 2 + 2f 2 F 2 + 29 2 sin 2 F j sin 2 6 sin 6d6df , (29) 



where M(R) is static mass functional (eq. (|TTD) and Inn , I^n > -^i^ are f° ur mo ~ 
ments of inertia of the skyrmion. In deriving equations (125|)-(|29D we have made use of the 
axial symmetry of the system. 

Defining now canonical conjugate variables in the body-fixed reference system 



T; 



and J-, 



dC 



one obtains the Hamiltonian of the system 

(I? + T 2 2 )4n 2) + (J! + JMi 2) + *{TJi + f 2 7 2 )7 ( 



H = M + 



M + 



.(12) .(12) 



y 1 ^ ) 



(12) ry 2 



47 



(33) 



T 2 7£? + J 2 e 2) + 2f • J7£) - 


^nn 






4 


7(i2) r (i2) r 




2" 





+ 



To 



47 



(33) 



u!2 



(30) 



(31) 



Here T is the isospin operator, J the spin operator. If in addition one considers states 
with K = 0, where K = T + J is the grand spin, and sandwiches the Hamiltonian between 
wave functions \T, T 3 ; 7, 7 3 = — T 3 ) in the body-fixed system of reference, the energy of 
the system takes the form 



E = M + 



7(12) , r(12) _ 2/ -(12)- 


(T(T + 


1) " If) 


4 


7(12) .(12) 


^j(12)\ 2 " 





+ 



T 2 

J 3 



47 



(33) 



(32) 



where we used the constraint T 3 = — 7 3 which is due to the axial symmetry of ansatz (|7|). 
The degeneracy on the quantum number 7 3 is partially lifted, i.e. the energy of the 
skyrmion only coincides for states with ±7 3 . For example, the following pairs of baryons 
have the same masses: 



>) 

Hi) 



\P, J 3 = - 
|A + , 7 3 = 
|A++,7 3 



-I) 



-I) 



and 
and 
) and 



I"-, 7 3 = 
|A°,7 3 
|A- J 3 



= 2 ) 
2/ 



(33) 



but the baryons in the states ii) and Hi) do not have the same masses. The above specified 
quantum numbers refer to the body-fixed system of reference. In contrast to a deformed 
skyrmion in an isotropic background (e.g. the vacuum or isotropic and isosymmetric 
nuclear matter) ||, our present case corresponds to a deformed nucleon or delta in an 
unisotropic background given by the nuclear density profile. Therefore, the projection onto 
good quantum numbers in the laboratory system of reference does not follow here from a 
simple rotation of the intrinsically deformed object. In fact, only the third components of 
angular momentum and isospin in the body-fixed system, 7 3 and T 3 , are good quantum 
numbers. 
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As a check, let us consider whether we can obtain the canonical spherically symmetric 
solutions from our deformed ones. First, it is easy to see that 



Km e 2) = lim ig? = lim iffi* = lim iff = \ , (34) 
0^# 0^# 0^# 0^9 

where I is the moment of inertia of the spherically symmetric hedgehog ansatz, 

oo 

/ = f sin 2 F (ffl + 4f 2p2 + 4 gin 2 F ^ d ~ _ (3^ 



Consequently 

r(12)r(12) _ fr(12)\ 2 / r(12) _ r(12)\ r(12) 

TP n 4- <?T^ ~ TP n r(12) .(12) ~ „ llm ^ ' W 

^ # ^ # ^ # 

and one obtains the well known energy expression for a baryon in the spherically sym- 
metric case 

E = M + T(T + 1) . (37) 

More generally, the following relation for all observables can be shown to hold: 

/ Expression in axially \ / Expression in spherically \ , . 

F e — > \ symmetric case J y symmetric case J 

This shows that our construction has the proper limit when all deformation parameters 
are set to zero. 



V. QUANTUM PROPERTIES OF SOLITONS 

A. Masses of quantized solitons 

Predictions for the masses of the quantized solitons are presented in Fig. [7], where 
we show the nucleon mass normalized to its free space value. Note that the distance R 
between the centers of the nucleon and the nucleus is normalized to the corresponding 
radius Ra of the nucleus for given baryon number A. We show the behavior of the nucleon 
mass only for one representative nucleus of the three classes, i.e. for 16 0, 56 Fe and 198 Au. 
The other nuclei from the same class exhibit a similar behavior. For the radii of the nuclei 
we use formula 
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/ R 2 p(R)d 3 R 

R A = (R 2 ) 1/2 = J — f — • (39) 

/ p(R)d 3 R 

The behavior of the nucleon mass in the nucleus is similar for all light, medium-heavy and 
heavy nuclei. When the density decreases, the mass of the nucleon increases approaching 
its free space value at the edge of the nucleus. We conclude that the mean effective 
mass of the nucleon in real nuclei will not decrease as much as ~ 40% in contrast to 
the case of homogeneous nuclear matter with normal density p = po- The finite nuclei 
approach for the medium essentially improves the results, especially for light nuclei, i.e. 
taking into account realistic mass distributions one obtains more reasonable values for the 
mass renormalization. For example, this renormalization is about ~ 17% for a nucleon 
placed in the center of 12 C or 16 O. In principle, one could now calculate an average 
medium modification in a given nucleus by simply folding the calculated effective mass 
with the nuclear density. We do not show the results of such calculations here because 
the calculation could be further refined. If one uses e.g. the shell structure of the nuclei, 
then all nucleons will be situated at some distances (within the shells) from their mass 
center and the mean effective value of the nucleon mass in such systems will be further 
reduced for nucleons in heavy nuclei. Further improvements of the results can be made 
by including a dilaton field |20|fjj2l||. We also note that a similar renormalization of the 
nucleon mass was found in Ref. 0. 

As it was noted in the previous section, the mass of the A states depends on the 
absolute value of the third component of isospin. We found small deviation from degen- 
eracy for delta states with J3 = —1/2 and J3 = —3/2^]. The difference of the masses 
AM A = M A +(J 3 = -1/2) -M A ++(J 3 = -3/2) of A in the states ii) and in) (see eq. (§3])) 
has a negative value near the center of the nucleus and after reaching its minimum begins 
to increase towards positive values reaching at some distance its maximal value and then 
smoothly drops to zero. For example, it has its minimal value AMa = — 0.28MeV at 
R = 0.84 fm and maximal value AM A = 0.916MeV at R = 2.447 fm for deltas in 12 C. For 
the case of 56 Fe, the corresponding values are AM&(R = 4.098 fm) = — 0.487 MeV and 
AM&(R = 7.806 fm) = 0.555 MeV. We remind the reader that all these numbers refer to 
states defined in the body-fixed system of reference as discussed before. 

For the relative mass difference M A — we obtained a behavior similar to the one 
of nucleon mass (cf. Fig. (]?])) and do not graphically present it here. The difference 
M A — Mjy decreases when the density increases and has its minimal value at the center of 
the nucleus. More specifically, (M A -M^)/(M A - M N ) has the value 0.867 for 12 C, 0.615 
for 56 Fe and 0.609 for 198 Au, in the respective center of the nucleus. Note that in free 
space M A — = 300 MeV in the current approach. As to the ratio of A and N masses, 
M A /M^, it is almost constant inside the nucleus and approximately equal to that ratio 
in free space. 



tt\y e remark that due to the neglect of the Fermi motion, these small effects should be consid- 
ered indicative. 
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B. Radii and intrinsic quadrupole moments of nucleons 



The calculation of the electromagnetic currents reduces to the calculation of the 
isoscalar and the isovector currents. In topological models and in particular in the Skyrme 
model, the isoscalar current is equal to the model-independent baryon number current 
(for a discussion on the role of vector mesons, see Ref. 



B<* = ^e^TrL v L a L p , (40) 

with e fJ,ua/3 the totally antisymmetric tensor in four dimensions. Similarly, the isovector 
current, which is the Noether current associated to the global transformations U(x) — > 
exp(iQ L )U(x) exp(—iQ R ) (Ql,r are 2x2 matrixes), has the form0 

^ = -^C^tt(L, + + -^Tr f{[L v , [Lp, L v ]\ + [R v , [R„ R„}}} , (41) 
lb loe 



where 



R, = Ud,U + ; C, 



1 , fi = 
a p , /i = 1,2,3 



Then, the operator of the electric charge of the nucleon is determined by the zeroth 
components of the isoscalar and the isovector currents, 

J e o m (r) = V i3) (r)f 3 + ±B (f), (42) 

(3) 

where Vq is the third component of the isovector current and T 3 is the third component 
of the isotopic spin operator. Sandwiching this current between corresponding nucleon 
states, we obtain the proton and neutron charge distributions. 

The normalized isoscalar root-mean-square radius along each axis, (r'^) I=0 , (normal- 
ized with the respect to the isoscalar charge) can be determined from the zeroth component 
of the baryon current and has the form 

OO 7T 

f f D a f 2 F f Q e sin 2 F sin QdOdr 

(J2\ -JLJ D fsm 2 fl, a = 1,2, 

^ ;/=0 ~ °?l ' ^ Q ~ 1 2cos 2 #, a = 3, ■ j 

2e 2 F 2 / F f e sin 2 Fsm Qd9dr 



o o 



where a (a = x, y, z) indicates the projection of the radius vector on the corresponding 
axis. Similarly, the isovector root-mean-square radii (normalized to the isovector charge), 
(r 2 ) J=1 , are determined by the zeroth component of the vector current Vq 3 ^ 



tt There is no summation over ju in the equation ( [4l|) while there is one over the index v 
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U*J it 

J J D a r 2 fsm9dBdr 

(r'*) I=1 = ° ° oo, , (44) 

2e 2 F2 J J fsmOdOdr 



o o 



where 

/ = (sinF) 2 [f 2 + 4F 2 + 4f 2 F| + 46 2 sin 2 F] . (45) 

One can see that the influence of the medium enters indirectly via the modified profile 
functions. Results for the in-medium radii are collected in Table |. We give the isoscalar 
(r 2 )]= and isovector (r 2 ) 1 /^ root-mean-square radii of the nucleons along the y and z 
axes inside the nuclei 16 O, 56 Fe and 198 Au, evaluated in the center of each nucleus and at 
those distances where the intrinsic quadrupole moment of the nucleons has its first and 
second extremum, correspondingly (see the discussion below). We note that due to the 
axial symmetry {r 2 )\^-i = (r 2 ) 1 /^.^. Clearly, these radii swell, as first pointed out in 
ref. |23[ in the analysis of data on quasi-elastic electron-nucleus scattering. However, this 
swelling is not isotropic and depends very sensitively on the distance from the center of 
the nucleus. The concept of a common swelling by one scale factor is certainly much too 
simple. Our results are therefore not in contradiction with the y-scaling analysis of ref. 
p4|| (see also the discussion in || on this issue). 

As noted before due to the deformation the nucleon can have a nonvanishing intrinsic 
quadrupole moment [^] 

Q int = ( J 3 \ J J cm (f) r 2 {3 cos 2 9-1) d 3 r| J 3 ) (46) 

that is characterized by the third component J3 of the angular momentum in the body- 
fixed frame of reference. The corresponding isoscalar and isovector intrinsic quadrupole 



moments follow from (f43|) and (|44| ) by appropriate projection 

OO 7T 

/ J (3L> 3 - 2)f 2 F f Q e sin 2 F sin QdOdr 



int (r>\ _ 



OO 7T ) 

2e 2 F 2 J J Ff Q e gin 2 p gin Q dedf 


(47) 



int /D^ _ 



J J(3D 3 - 2)f 2 f sin OdBdf 



Qfli(R) 



2e 2 F 2 J J f s i n Qd6dr 







Consequently, the proton and neutron intrinsic quadrupole moments follow from 

= \(Q?1 + QTli) , (ff = liQflo - QTli) ■ (48) 
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The dependence of the isoscalar and isovector intrinsic quadrupole moments as a function 
of R/Ra is shown in Fig. ||. Again we present results only for the nuclei 16 0, 56 Fe 
and 198 Ag. For the other nuclei the dependences are similar. One can see that intrinsic 
quadrupole moments are small. This smallness of the intrinsic quadrupole moments means 
that nucleons in the nucleus are weakly deformed. The intrinsic quadrupole moments has 
two extrema. Near the center of nucleus they have negative values and a corresponding 
minimum at that distances where density changes are ~20% for 16 and ~7% for 56 Fe 
and 198 Au, respectively. Then with increasing density the intrinsic quadrupole moments 
become positive and have a maximum at the distance where the density of the nucleus has 
fallen by ~ 65 4- 70%. Note also that the isoscalar intrinsic quadrupole moment is about 
two times smaller than the isovector one. Consequently near the center of the nucleus 
protons are flattened along the axis z (oblate) and and near the border of the nucleus they 
are extended (prolate). In contrast neutrons near the center of nucleus slightly prolate 
and near the border of nucleus oblate. Consequently, at some intermediate distance the 
nucleons must have spherical form. 

Finally, let us recall that the observable static quadrupole moment in the laboratory 
system of an intrinsically deformed baryon in an isotropic background is linked to its 
intrinsic quadrupole moment by the Wigner-Eckhardt theorem as follows 

Q lab = (J J 3 2 | JJ 3 ){J ± J 2 | J ± J) Q int 

3j|_j (J+ l) JgJ-1) int 

" (J+l)(2J + 3p (J+l)(2J + 3)^ ' { yj 

where the relations are valid for maximally projected states with |J 3 | = J. Thus for a 
deformed baryon of total angular momentum J = 1/2 the quadrupole moment in the 
laboratory system is zero for its ground state, i.e. for the nucleon. It would be interesting 
to investigate the transition quadrupole moments to excited states, extending the work 
of ||. As mentioned before, our case corresponds to a deformed baryon in an unisotropic 
background. A simple projection to the laboratory system is therefore not possible. 



VI. CONCLUSION 



In summary, we have considered properties of a nucleon embedded in finite nuclei 
making use of the Skyrme model. Deformation effects are taken into account by the 
distortion of chiral profile function under the action of the nuclear medium, parameterized 
in terms of the functionals Eq. flTT|) . Nucleons and deltas emerge as quantized solitons after 
performing an adiabatic rotation in co-ordinate and isospace which takes into account the 
axial symmetry of the system (leading to two different moments of inertia, I z ^ I x = I y ). 

The main results of this study can be summarized as follows: 

(i) The density dependence of the nucleon mass shows a more realistic behavior than in 
the case of a uniform density as e.g. for homogeneous nuclear matter. The effective 



§§The first Clebsch-Gordon coefficient refers to the body- fixed frame, while the second one 
refers to the laboratory frame. 
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nucleoli mass has its minimum at the center of the nucleus and approaches its free 
space value at the surface of the nucleus. 

(ii) The nucleons in finite nuclei are weakly deformed. They acquire a small intrinsic 
quadrupole moment which, however, is strongly dependent on the distance from the 
center of the nucleus, e.g. the proton (neutron) deformation changes from a oblate 
(prolate) to a prolate (oblate) form as it is moved toward the surface of the nucleus. 

(iii) Similarly, there is a direction-dependent swelling of the isoscalar and isovector root- 
mean-square radii. We have stressed that the concept of a uniform swelling factor 
is too simple a concept to apply to real nuclei. 

(iv) As consequence of the axial symmetry of the system, the | J^\ = 1/2 A states (A°' + ) 
and the | J 3 | = 3/2 ones (A~' ++ ) have slightly different masses in finite nuclei. 

Let us briefly give an outlook for further studies based on these results. As it was 
mentioned above this approach can be used for model investigations of the shell structure 
of nuclei and to study the properties of a remote nucleon in such a system (halo nuclei). In 
the latter case the nucleon in the finite nucleus should not be considered as one within the 
nucleus but weakly coupled to the nucleus at large distances. Also, the model naturally 
allows to study nucleon properties in the lightest nuclei, especially the electromagnetic 
structure of nucleons in 3 He, which is of current experimental interest. Work along these 
lines is under way. 
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TABLES 



Element (rg)j^ [fin] (r^g, [fm] (rg)^ [fin] (r z 2 )}g [fm] fl [fm] p(fl)/p(0) 

Center of the nucleus 



16Q 


0.437 


0.437 


0.699 


0.699 





1 


56 Fe 


0.520 


0.520 


0.811 


0.811 





1 


198 Au 


0.523 


0.523 


0.815 


0.815 





1 








First extremum 








16Q 


0.427 


0.425 


0.684 


0.682 


0.74 


1.20 


56 Fe 


0.499 


0.493 


0.784 


0.778 


2.97 


0.93 


198 Au 


0.503 


0.495 


0.788 


0.780 


5.45 


0.93 








Second extremum 






16Q 


0.404 


0.411 


0.651 


0.661 


2.13 


0.35 


56 Fe 


0.433 


0.442 


0.693 


0.708 


5.02 


0.32 


198 Au 


0.435 


0.444 


0.696 


0.711 


7.46 


0.31 








Free space 









0.391 0.391 0.630 0.630 



TABLE I. The isoscalar (r 2 ) I=0 and isovector {r 2 ) I=1 nucleon root-mean-square radii along 
the y and z axes inside various nuclei. The values of these quantities are given for the center of 
the nuclei and at those distances where the nucleon intrinsic quadrupole moments have their first 
and second extremum values, correspondingly. R is the distance between centers of the nucleus 
and skyrmion topological center and the last column is the ratio of the density of nucleus p(R) 
to its value at the center p(0). We note that (r 2 ) 1 /^^ = (^)/= -i- ^or comparison we also give 
the free space values of the root-mean-square radii. 
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FIGURES 




FIG. 1. Location of a skyrmion inside the given nucleus. O is the center of the nucleus 
and O' is the topological center of the skyrmion. We define by \R\ the distance between these 
centers. The direction of R coincides with directions of axes z and z 1 . 
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2.0 2.5 3.0 

R.fm 



FIG. 2. Mass dependence of a skyrmion embedded in a given nucleus (ratio of soliton mass 
in the nucleus to its free space mass, Mg/Ms) as a function of the distance R between centers of 
nucleus and skyrmion (band 1), the ratio of the density at distance R to its value in center of the 
nucleus, p/p(0) (band 2) and the dependences of the deformation parameters 5rs/rs (band 3), 
71 (band 4), 72 (band 5), 5\ (band 6) on distance R. Note that 5\ is multiplied by factor 10. 
Solid lines represent these dependencies for a skyrmion in 12 C and the dashed ones for that one 
in 16 0, respectively. We note that the ordinate scale of the top and bottom figures are different. 
For definitions of 8rs/rs, 71, 72 and 6\ see eqs. (HD, (ID and ©• 
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1.10 




-0.255 hin" 

0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 

R.fm 



FIG. 3. As in the Fig. § but for a skyr mion embedded in Ca (solid lines) and 56 Fe (dashed 
lines), respectively. For notations see Fig. & 
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FIG. 4. As in the Fig. || but for a skyrmion embedded in 198 Au (solid lines) and 208 Pb 
(dashed lines), respectively. For notations see Fig. |2[ 
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a) 




-3-2-10 1 2 3 

Z 



b) 

FIG. 5. Baryon charge distribution Bq of the skyrmion embedded in 12 C in the yz plane 
(a) and its projection (equal baryon charge distribution areas) onto this plane (b), where y, z 
are dimensionless variables. The distance from the center of the nucleus R is equal to 1.266 fm. 
This corresponds to maximum value of the leading deformation parameter 71. The center of 
the nucleus is situated on the negative side of the axis z. The strength of the baryon charge 
distribution corresponding to the different area is represented by light /dark shading. The baryon 
density is lower if the area is darker. The values of the deformation parameters at this place are 
§r s = -0.0476 fm, 71 = -0.0763, 72 = 0.028, 73 = -0.0039, Si = 0.0029 and 6 2 = -0.0003. 
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a) 
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FIG. 7. Mass dependence of a nucleon embedded in various nuclei as a function of the 
fraction R/Ra, where R is the distance between the centers of nucleus and nucleon and Ra 
is the radius of the nucleus (see eq. (p9|)). The solid line represents this dependence for the 
nucleon in 16 0, while the dashed and dotted lines represent this dependence for 56 Fe and 198 Au, 
respectively. 
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FIG. 8. The R/Ra dependence of isoscalar (top figure) and isovector (bottom figure) 
intrinsic quadrupole moments of a nucleon embedded in various nuclei. For the notations see 
Fig. ft 
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